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Abstract 

We study properties of compact stars with the deconfinement phase transition 
in their interiors. The equation of state of cold baryon-rich matter is constructed 
by combining a relativistic mean-field model for the hadronic phase and the MIT 
Bag model for the deconfined phase. In a narrow parameter range two sequences 
of compact stars (twin stars), which differ by the size of the quark core, have been 
found. We demonstrate the possibility of a rapid transition between the twin stars 
with the energy release of about 10 52 ergs. This transition should be accompanied 
by the prompt neutrino burst and the delayed gamma-ray burst. 

PACS: 14.65.-q, 26.60,+c, 97.10.-q, 98.70.Rz 

1 Introduction 

Our main goal in this paper is to study properties of compact stars composed of strongly 
interacting matter undergoing the deconfinement phase transition. We construct the 
equation of state (EoS) at finite baryon density ps by combining two popular models of 
the baryon-rich matter. For the hadronic phase we use the NLZ version of the relativistic 
mean-field model |I[ which gives very good description of the saturation properties of cold 
nuclear matter. The deconfined phase is described within a simplified version of the MIT 
Bag model ||. The possibility of a phase transition between these two phases has been 
carefully studied by applying the Gibbs conditions for charge-neutral and /3-equilibrated 
matter. Indeed, we have found a first order phase transition and determined characteris- 
tics of the mixed phase. 

Properties of hybrid stars containing both the hadrons and quarks have been already 
studied by using a large number of models (see e.g. Refs. [§, [5|, §, [7j). One may wonder, 
what is new in our work? First of all, we want to note that this paper presents only a 



small part of our more comprehensive study || where we systematically analyze many 
other models of hadronic and quark phases. It is interesting that in most cases we do 
not find any phase transition between the two phases. The combination of the NLZ and 
MIT Bag models represents one of a few exceptional cases, when the deconfinement phase 
transition is predicted in stellar interiors. Moreover, we find two families of compact stars, 
twin stars || [K| |ll|] , which differ by the size of the quark core. This opens the possibility 
of a catastrophic rearrangement of the twin star from one to the other configuration with 
a release of energy of about 10 52 erg. In this paper we present these results and discuss 
their observational consequences. 

2 Properties of matter in compact stars 
2.1 Hadronic phase 

There is no doubt that hadrons, mesons and baryons, are correct degrees of freedom for 
modeling strongly-interacting matter at low densities. From nuclear phenomenology we 
know that atomic nuclei can be well described in terms of interacting nucleons. Therefore, 
we believe that the hadronic phase should be stable at least up to the saturation density 
of nuclear matter po — 0.15 fm -3 . Field-theoretical models, where nucleons interact with 
mean meson fields, are proved to be very successful in describing saturation properties of 
nuclear matter as well as properties of finite nuclei. Here, to calculate the EoS of hadronic 
matter we use a non-linear version of the relativistic mean-field model known as the NLZ 
model Compared wih the original version it is generalized by including hyperons (Y) 
and hyperon-hyperon (YY) interactions, as proposed in Ref. ||12|| . We call this model 
NLZY. 

The original Lagrangian density for the NLZ model (without the YY interaction) can 
be written as [[|] (here and below we use the units % — c — 1) 

£ = J2^b (ip-m B ) 1>b + ^ad^a - ^rn 2 a a 2 - V - V - ^u^u^ + ^m* w"^ 

~ JP^Ptxv + \™ffiPv + (daBCT + guBU^ + gpP^Ts) ^B , (1) 

4 2 B 

where the sum runs over all the baryons B=p, n, A, H 0,± , H 0, ~. In the above Lagrangian 
a, u and p are the iso-scalar scalar a, the iso-scalar vector u and the isovector vector p 
meson fields respectively. In Eq. (|l]) and p^ y denote, respectively, the field strength 
tensors for the uo and p meson fields. 

Originally this model was designed for the nucleonic sector and it failed to reproduce 
the observed strong AA attractive interaction. This defect can be removed by adding 
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two new meson fields with hidden strangeness, namely, the iso-scalar scalar a* and the 
iso- vector vector 0, which couple to hyperons only |I2"| |. These fields can be identified 



with the Jo (975) and 0(1020) mesons. The corresponding Lagrangian is given by 
C YY = ~ (d»a*d»a* - mla* 2 ) - + \m\^^ + E (9«*yct* + g^ly) 

where index Y runs over hyperons only. The mean meson fields are found from the 
Euler-Lagrange equations. 

The nucleon coupling constants are chosen from the fit of the finite nuclei properties. 
The vector coupling constants of the hyperons are chosen according to the SU(6) sym- 
metry and the hyperonic scalar coupling constants are chosen to reproduce the measured 
values of the corresponding optical potentials. Below we use the set of model parameters 
suggested in Refs. 0, |12| - 

Within the mean-field approximation, the pressure and energy density of static and 
homogeneous baryonic matter can be easily calculated from the above Lagrangian: 

H 1 2 2 , ^ 3 , C 4 , ^ 2 *2,1 22, 

e = 2 m - a + 3 + 4 a +2 m<T * a +2 m ^o+ 

+ \m 2 pP l, Q + \m\<& + E ^ f F dkk 2 J¥^tf , (2) 

n H ^ 2 2 ^3 C 4 1 2 *2 . 1 2 2 1 

P = — m„a a a m„,a H — m, ,uj n + 

2 ° 3 4 2 T 2 w0 

loo lo,o x — f k F k 4 



+ \rny ofi + \m\<& + E^/ " dk > (3) 



where m* B = rris — g a B& — 9a* b&* is the effective mass, vb is the degeneracy factor and 



kp = y fi B — m* B is the Fermi momentum of the baryon species B. 

In order to have a complete description of the /3-equilibrated matter one should also 
include the leptons. In both the hadronic and quark phases their contributions to energy 
density and pressure are given by the well known formulae of ideal Fermi gas. 



2.2 Deconfined phase 

As follows from a simple geometrical consideration, nucleons begin to overlap at densities 
Pb ~ {^ Kr N /3) _1 — 3po for the nucleon radius — 0.8 fm. Such densities are surely 
reached in the interiors of compact stars. Of course, this argument does not tell anything 
about the character of transition from hadronic to quark-gluon degrees of freedom. Below 
we follow the common practice of using two different models for these two phases. Namely, 
the deconfined phase is described within a simple version of the MIT Bag model 0], 
considering it as a mixture of free Fermi gases of u, d, s quarks in a bag with an additional 
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energy density B (the bag constant). Within this model the energy density and pressure 
of cold deconfmed matter are written as 

J2 ^ /*' dkk?y/m} + k* + B, (4) 

_ , Z7T Z JO v 



f=u,d,s 



P Q = E r'dk , = -B, (5) 



where k F = Jfj/j — mj is the Fermi momentum of quarks with flavor /. For each flavor 
we choose the degeneracy factor Uf = 2 (spin) x 3 (color) = 6 and take the following values 
of quark masses: m u = 5 MeV, = 10 MeV and m s = 150 MeV. 

2.3 Conditions of local equilibrium 

In /3-equilibrium, the chemical potential of any particle species i can be expressed as 

Hi = kHb + QiHe , (6) 

where hi is the baryon number of the species i , denotes its charge in units of the 
electron charge, fj,b an d He are the baryonic and electric chemical potentials, respectively. 
Here and below we assume that neutrinos can freely escape from the star. Eq. (§) means 
that only reactions conserving charge and baryon number are allowed. On the other 
hand, strangeness is not conserved because strangeness-changing reactions are generally 
much faster than a characteristic time of the star evolution. Two independent chemical 
potentials, /ij, and found by fixing the baryon and electric charge densities: 

Pb = E h *Pi ' Pe = E ftPi > ( 7 ) 

i i 

where pi is the number density of the particle species i . It is obvious that stars must be 
electrically neutral on a macroscopic scale, i.e. p e = 0. 

We assume that deconfinement is a first order phase transition which, in general, 
should produce a mixed phase (MP) between the pure hadronic phase (HP) and pure 
quark phase (QP). At zero temperature the MP should follow the Gibbs conditions 

P H (H b ,He) = P Q (Hb,He), (8) 
Hb = Hb = H? » ( 9 ) 

A*e = A*f = H Q e- (10) 

According to Eq. the baryon chemical potential \ib equals the neutron chemical 
potential \i n and p e is equal to the electron chemical potential. At given //& and /i e , 
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the quark chemical potentials are found by using the formulae /i u = (fib — 2/i e )/3 and 

The volume averaged energy density in the MP can be written as 

e = (l-\)e H (fx b ^ e ) + \e Q (fx b ^ e ), (11) 

where A = Vq/V is the volume fraction of quark phase. In the case of two chemical poten- 
tials one can only construct the MP by adopting a generalized (global) charge neutrality 
condition [13|], when the net positive charge of one phase is compensated by the negative 
charge of the other phase: 

(l-A)pf(^,/x e ) + Ap«( yUb , yUe ) =0. (12) 
This condition is assumed in our calculations of hybrid stars presented below. 




Figure 1: Pressure as a function of fi b and /i e for the hadronic phase (wired surface) as 
predicted by the NLZY model and for the quark phase (grey surface) as predicted by the 
MIT Bag model (5 1 / 4 = 180 MeV). The thick line represents the charge neutral matter. 
Its portion between points C and D corresponds to the mixed phase. 
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2.4 EoS and composition of matter 

Figure [l] shows the pressure surfaces for the pure HP and pure QP (for bag constant 
B l / A = 180 MeV) as functions of chemical potentials [L\, and p e . It is important that 
these two surfaces intersect, which is not the case for most of the other models ||. In 
this case one can construct a MP connecting the pure phases along the intersection line. 
The thick curve on the surfaces shows pressure of charge neutral matter in /3-equilibrium. 
The lower part of this curve (from low pressure to point C) corresponds to the pure HP. 
The MP starts at point C corresponding to baryonic density ps — 1.5 po , and ends at 
point D corresponding to ps — 5.1 po • At higher densities the matter is composed purely 
of quarks. Our calculations show that at B 1 ^ > 190 MeV the pressure surfaces of the 
HP and QP do not intersect at all and the MP can not be defined by the Gibbs rules. 

The particle composition in these three regions is presented in Fig. [| It is interesting 
to note that the only hyperon species which survive in the MP is A-particle which appears 
at densities 2.8p < Pb < 5.1p . All the negatively charged hyperons are suppressed as a 
result of the charge neutrality. The strangeness content of matter rapidly grows in the MP 
and reaches the asymptotic value of 1/3 in the QP. On the other hand, leptons (mainly 
electrons) are practically extinguished at ps>3po- 




Figure 2: Relative particle abundances versus baryon density p^ as predicted by the 
combined NLZY and MIT Bag (E 1/4 = 180 MeV) models. The region 1.5 < p B /po < 5.1 
corresponds to the mixed phase. 
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3 Properties of compact stars 



Below we study properties of static spherically symmetric stars. Under assumption that 
the matter may be treated as an ideal fluid, the star structure can be found by solving 
the TOV equations [14]. For a given EoS, P = P(e) , and a fixed central baryon density 
Pc = Pb{t = 0) we integrate the TOV equations from the center of the star up to its 
surface r = R. The star radius R is determined from the condition P(R) = 0. The details 
of numerical calculation as well as the online program can be found on the internet home 



page [HJ] of one of the authors (M.H. 




Figure 3: The mass-radius relations for different families of stars. Open dots correspond 
to the maximum mass. The beginning and the end of the mixed phase in hybrid stars are 
shown by full squares and dots. 

In Fig. ^ we present the mass-radius relations calculated for several values of B as well 
as for pure hadronic and quark stars. By open dots we show the critical configurations 
with highest possible masses. The pure hadronic curve are calculated assuming the NLZY 
EoS at all densities. The part of this curve corresponding to M > 1.6M is omitted due 
to appearance of densities with negative effective masses of baryons [P^j . The beginning 



of the MP at the star center is marked by full squares whereas its end and, therefore, the 
beginning of the pure QP, is marked by full dots. As seen in Fig. ^|, the model predictions 
depend strongly on the bag constant. Stable hybrid stars are possible only for low bag 
constants, B 1 ^ < 180 MeV. For higher B stars become unstable even before they reach 
high enough density for the formation of a pure deconfined phase in their interiors. On 



7 



the other hand, pure quark stars are possible at any B . 

A very interesting feature is found in a narrow interval of B around 5 1/4 = 180 MeV, 
where a second sequence of stable hybrid stars appears. Their properties can be summa- 
rized as follows: 

• The first sequence of stars can only reach maximum central baryon densities of 
about 0.77 fm -3 which is just at the end of the MP. As one can see from Fig. || 
stars of this sequence have masses below M max = 1.36 M and radii R > 13 km. 
These stars are mainly composed of hadronic and MP matter with admixture of the 
A-hyperons. The calculation shows that stars near the maximum mass contain a 
tiny core of pure quark phase. 

• Stars of the second sequence have considerably higher central densities, 0.9 fm~ 3 < 
p c < 1.53 fm -3 . Their masses lies in a narrow interval, 1.35 < M/M Q < 1.38 and 
their radii are noticeably smaller, 10.2 km < R < 12.3 km. These stars are mainly 
composed of quark and MP matter, surrounded by a rather thin layer of HP and 
the nuclear crust. 

• The two sequences of compact stars are separated by unstable region corresponding 
to the interval of central densities 0.78 fm -3 < p c < 0.9 fm -3 . 

From our more comprehensive study |8j involving many other models of EoS we con- 
clude that two sequences of compact stars, sometimes called twin stars, appear only in 
very exceptional cases. There is no guarantee that this picture is the one that corre- 
sponds to the reality. But we find it interesting to study its possible consequences for the 
dynamics and observable signatures of compact stars. 

4 The twin star collapse 

One can imagine the situation when a compact star on the first sequence has a mass close 
to M max . If this star is a member of binary system, then its mass and, therefore, its 
central density may grow due to accretion from a companion star. Eventually the mass 
will exceed the maximum value when the star becomes unstable with respect to radial 
compression. Usually it is assumed that this loss of stability leads to the collapse into a 
black hole. However, our calculations open another possibility: the collapse into the twin 
star on the second sequence. 

Let us assume that no matter is ejected during this process, i.e. the total baryon num- 
ber, Nb is conserved. Figure [| shows the gravitational mass of the star as a function of its 
total baryon number. The star from the first sequence which reaches the maximum mass 
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Figure 4: Mass as a function of baryon number for different families of stars. The inset 
shows the enhanced region of twin stars. The twin star collapse corresponds to the 
transition from A to B. 

(point A) will collapse to its twin star. The latter is the corresponding star on the second 
sequence, i.e. the one which has the same total baryon number (point B). The difference 
in energy between these two stars, AE, is given by the difference in their gravitational 
masses. In the case considered here, the released energy is AE ~ 6 x 1O~ 3 M ~ 10 52 erg. 
This amount of energy should finally be emitted in one or the other way. 

The baryonic density profiles of the twin stars A and B are compared in Fig. ||. One can 
see that an extended MP is present in both cases but the star B has much larger quark core 
as compared to the star A, where this core is only marginally present. So, the instability 
develops practically at the end of the MP region. In other words, this shows that stars 
with small quark cores are unstable. In the considered example the new equilibrium state 
appears only when the core radius exceeds about 4 km. This reminds the well known 
result from the theory of nonrelativistic stars with a density jump inside. Namely, if pi 
and pi are baryon densities just below and above this jump, then small dense core becomes 
unstable if p\j p2 is larger than a certain critical value (3/2 for incompressible matter 
Of course, in our case there appears no jumps of density, although its gradients are large 
near the core boundary. 

This kind of instability associated with a first order phase transition in hadronic matter 



was first studied in Ref. [171 and later on in Refs. 18, 191. It was shown that the transition 
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Figure 5: Baryon density profiles of the twin stars A and B as defined in Fig 4. 

from a small-core to a large-core configuration proceeds via large-scale damped oscillations 
around the new equilibrium state. We expect that the collapse of twin stars may proceed 
in a similar way. After reaching the critical state A the star looses stability and enters the 
stage of catastrophic rearrangement. Since the QP in the core has higher density than 
the replaced MP, surrounding layers of the star will acquire collective inward motion. 
Due to inertial effects, the star will overshoot the new equilibrated state B and rebound. 
This will give rise to the damped oscillations around this state. As simple dimensional 



estimates show [[17]], the initial implosion and following oscillations are characterized by 
the millisecond time scale. 

Such star rearrangement may involve many interesting processes. First of all, large 
amount of hadronic matter will be transformed into deconfined phase. Since strangeness 
and leptonic contents are very different in two phases, this transformation will require 
not only strong-interaction reactions but also weak processes changing strangeness and 
lepton numbers. Such reactions will inevitably produce large numbers of neutrinos. For 
instance, to achieve chemical equilibrium in the deconfined phase one needs reactions 
e + u — > d + v and u + d — > s + u which lead, respectively, to emission of neutrinos and 
production of additional strange quarks. We believe that the time scales of these weak 
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processes are much shorter than the characteristic times of the star collapse. 

The total number of neutrinos, N v , produced during the transition process can be 
estimated as follows. From Fig. [5] one can see that the final star has the quark core 
with baryon density p core ~ 5po and radius i? C ore — 4 km. Thus, the core contains 
N core ~ 2 x 10 56 baryons. They were initially in the HP at lower density p# ~ 2po . 
According to Fig. || the fraction of electrons in this phase was about 10% and practically 
no electrons were in the QP. Therefore, the difference between the number of electrons 
in configurations A and B is AN e ~ 2 x 10 55 (the direct numerical calculation gives 
A7V e = 1.6 x 10 55 ). Because of lepton number conservation, all these electrons should 
be transformed into neutrinos and therefore, N v = AN e ~ 2 x 10 55 . This number is 
about 1% of the total baryon number of the star N B ~ 1.7 x 10 57 . The neutrino energies 
will cover a broad range up to a maximum value E™ ax = p e + p n — p^ . One can estimate 
£jma.x assum i n g that initially the deconfined phase was produced in the nonequilibrium 
state with the flavor composition corresponding to the HP. The latter consisted mainly 
of neutrons, with the quark structure udd , as well as about 10% admixture of protons 
and electrons. Taking p d ~ 2p„ ~ 10p and p e ~ 0.5p we get E™* ~ 150 MeV which 
corresponds to the mean neutrino energy ~ 100 MeV. Multiplying this energy by N v one 
obtains that a significant fraction ( > 30%) of the released energy will be carried away by 
the prompt neutrino burst. 

We expect that the remaining energy will be transformed into heat. Nonequilibrium 
processes during the phase transformation as well as viscosity effects might be responsible 
for dissipation of collective kinetic energy and eventually for damping of oscillations. 
Assuming that thermal energy AE? is initially dissipated in the quark core, one can 
estimate its temperature T from the relativistic Fermi gas formula 

7r 2 T n 2 AE T . . 

*± Pg ^'core 

where p g = Pb/3 ~ 0.4 GeV is typical quark chemical potential at p B ~ p core . For 
AEt ~ 7 x 10 51 erg Eq. ([13]) gives To ~ 40 MeV. At later times the heat wave will prop- 
agate through the whole star and finally will produce photons and, possibly, electron- 
positron pairs at its surface. The emission temperature, T em , can be estimated by as- 
suming that the thermal energy is distributed over the whole stellar matter, which is 
nonrelativistic outside the core. Then one obtains T em ~ 2AE T /3N B ~ 2 MeV. These 
predictions give us a ground to think that the discussed mechanism may serve as the en- 
gine for Gamma- Ray Bursts (GRB) p0| . Another obvious prediction is that the discussed 
rearrangement of the star will lead to a significant change in its moment of inertia. In a 
rotating star this will result in a super-glitch phenomenon [PTf] . 
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5 Conclusions and outlook 



It is shown that within a realistic hybrid (NLZY-MIT Bag) model two sequences of com- 
pact stars (twin stars) are possible. We demonstrate that their interiors differ mainly by 
the size of pure quark core. The energy difference between two twin stars with the same 
baryon number is about 10 52 erg which is approximately 1% of their total energy. The 
transition between twin stars may be triggered by accretion of mass from a companion 
star or by some other processes leading to increase of the central density above the critical 
value. After that a catastrophic rearrangement of the star begins. It will first collapse 
and then oscillate around a new equilibrium state. 

Our estimates show that in the course of this process the quark core can be heated up 
to temperature of about 40 MeV which is comparable to the one in supernova explosions. 
Our main prediction is that the transition between twin stars will produce a prompt 
burst of neutrinos (with energies of about 100 MeV) followed by a gamma ray burst 
(with photon energies of about 1 MeV). These features of the twin star collapse make it 
a potential candidate for the GRBs. On the other hand, there are some characteristics 
of GRBs, e.g. beaming of radiation, which is impossible to explain without additional 
assumptions invoking rotation and strong magnetic fields. This problem needs further 
study. 

It is interesting to note that the presented model predicts maximum masses of compact 
stars, (1.35 — 1.38)M Q , which agree well with observed masses of pulsars [|[]. Moreover, 
the twin stars of the second sequence have nearly the same masses in the broad range of p c . 
If pulsars have quark cores, this may explain why their measured masses are clustered 
near the mass of about 1.4M . 

Other phase transitions may take place in the interiors of compact stars, e.g. kaon 
condensation, color superconductivity etc. In Ref. || it was shown that formation of 
metastable hyperonic matter can also lead (in a certain range of parameters) to the twin 
star solutions. The comparison of twin star properties predicted by different models will 
be given in a separate work ||. 
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